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Abstract — We consider a linear feeder connecting multiple 
distributed loads and generators to the sub-station. Voltage is 
controlled directly at the sub-station, however, voltage down 
the line shifts up or down, in particular depending on if 
the feeder operates in the power export regime or power 
import regime. Starting from this finite element description 
of the feeder, assuming that the consumption/generation is 
distributed heterogeneously along the feeder, and following the 
asymptotic homogenization approach, we derive simple low- 
parametric ODE model of the feeder. We also explain how 
the homogeneous ODE modeling is generalized to account for 
other distributed effects, e.g. for inverter based and voltage 
dependent control of reactive power. The resulting system of 
the DistFlow-ODEs, relating homogenized voltage to flows of 
real and reactive power along the lines, admits computationally 
efficient analysis in terms of the minimal number of the feeder 
line "media" parameters, such as the ratio of the inductance- 
to-resistance densities. Exploring the space of the media and 
control parameters allows us to test and juxtapose different 
measures of the system performance, in particular expressed in 
terms of the voltage drop along the feeder, power import/export 
from the feeder line as the whole, power losses within the 
feeder, and critical (with respect to possible voltage collapse) 
length of the feeder. Our most surprising funding relates to 
performance of a feeder rich on Photo Voltaic (PV) systems 
during a sunny day. We observe that if the feeder is sufficiently 
long the DistFlow-ODEs may have multiple stable solutions. The 
multiplicity may mean troubles for successful recovery of the 
feeder after a very short, few periods long, fault at the head of 
the line. 

I. Introduction 

Motivation. Typical feeder of a distribution system forms 
a line bringing power to a large number of sequentially- 
connected customers. Electric characteristics of the feeder 
are normally uniform while a snapshot of consumption varies 
from neighbor to neighbor even if these are statistically 
similar customers. We are interested to describe power 
flows and voltage profile along the line, where the latter 
is maintained constant only at the head of the line. When 
feeder consists of standard consumers, which extract both 
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real and reactive power from the system, the voltage sages 
down the line showing the largest deviation from the nominal 
value at the end of the line. Standard design principe for 
this type of long but simple line is to have it limited in 
length, so that the typical voltage drop at the end of the 
line would not exceed the 5% which can still be tolerated. 
However, the consumption pattern vary thus leading from 
time to time to an incidental off-bounds drop or raise of 
voltage. One may also discuss lines of varying length, asking, 
in particular, what is the maximal possible length of a 
feeder which supports (within given model of consumption) 
a valid power flow solution? Moreover, one is interested 
to generalize, thus discussing not only feeders consuming 
electricity but also producing, i.e. these built of distributed 
generators, e.g. representing small-scale PV-systems. Such 
PV-enforced customers are consumers when it is cloudy 
turning into producers on a sunny day [1]. Line consisting of 
distributed producers may inject power in the transmission 
system (if real losses along the line are small), also leading 
to voltage increase at the end of the line with respect to the 
level set at the head of the line. 

Related Work. Discussion of the static Power Flow (PF) 
equations in the case of a distribution feeder was pioneered 
in [2], [3]. It was shown, that it is natural and convenient to 
state the PF equations in this one-dimensional case in terms 
of voltages, and real and reactive powers flowing through 
segments of the line. In this study we will also adopt the so- 
called DistFlow formulation of the PF equations from [2], 
[3]. The main technical thread of the present manuscript is 
related to the DistFlow equations considered in the limit of 
large (continuously many) consumers forming the feeder. In 
particular we derive and analyze the Distflow ODE version of 
the original DistFlow. ODE of this type was first introduced 
and analyzed in [4], where the main focus was on recovering 
the classic power engineering phenomena of voltage collapse 
[5], [6], [7], [8], [9], [10] in the continuous setting of a long 
statistically homogeneous feeder. Note also that a related 
description of the power system as a media was made in 
[11], [12], where propagation of electro-mechanical waves 
through transmission system was introduced and discussed. 
(This type of waves do not develop in the transmission 
system, simply because of the lack of big inertial machines in 
the distribution. Also, and for the same reason, the electro- 
mechanical waves are not generated within the feeder.) In 
addition to following the Distflow ODE thread of [4], this 
paper also continues and extends analysis and approach of 
[13], [14], [15] in what concerns inverter-based distributed 
voltage control in the distribution feeder. 



Our Contribution. The task of this manuscript is two 
fold. First, and taking the algorithmic/computational stance, 
we aim at developing low-parametric (spatially homoge- 
nized) ODE modeling of a feeder consisting of many similar 
consumers. In some of the most general regimes, e.g. in the 
regime correspondent to constant consumption/production 
of power along the feeder, ODE model posses symmetry 
which allows funding power flows and voltage profile most 
efficiently. On the other hand, and once the homogenized 
ODE model is stated and justified, we aim to analyze power 
flows in the long distributed feeder, paying a special attention 
to discovering qualitative phenomena of interest for operation 
and control of the feeder in the interesting regimes. To 
complete the second task we analyzed distributed feeder in 
different regimes. For example, we study the feeder in the 
regime of a standard distributed consumption of real power 
and analyze different voltage control schemes associated with 
the freedom in choosing consumption or injection of reactive 
power through the distributed inverters. We also analyze 
future possibility of using feeder rich on PV systems in the 
regime where the feeder would provide export of power to 
the transmission (high voltage) part of the grid. Our most 
surprising conclusions apply to this case. We have found that 
this export regime may be extremely dangerous as generating 
more than one stable solutions. Then the healthy (normal 
voltage) coexists with possibly many low-voltage branches. 
Normally the healthy branch is to be realized, however a 
significant perturbation, e.g. a short fault at the head of the 
line, may force the system to transition into a low-voltage 
branch where the system would stay till the inverters are 
disconnected/tripped (or burned). The dangerous states are 
complex, in particular showing oscillations of power flow 
directions along the feeder. We also find that some previously 
suggested reactive control, even though improving the quality 
of the healthy solution, does not remove the dangerous 
degeneracy. 



The material in the rest of the manuscript is organized 
as follows. We state the standard DistFlow formulation 
and show how the homogenization approach results in the 
DistFlow ODEs in Section |TTJ We review recently suggested 
reactive flow/voltage control schemes in Section III Section 



IV is devoted to description of how the ODE formulation 
of the power flow equations allows re-scaling and even 
further reduction in the number of key characteristics and 
consequently efficient computations of the Cauchi (initial 
value) type. We present simulations of the feeder working in 
the consumption/generation, control/no-control regimes and 
discuss the results in Section [V] Section [VI] summarizes the 
material and suggests a path forward. 



II. Homogenization: from Finite-Element 
DistFlow to DistFlow-ODEs 

The DistFlow equations of Baran, Wu [2], [3] have the 
following form 



Pk+l - Pk = Pk - rk 



P k+Qt 



Qk+i-Qk = qk-x k 



Pk+K 



v 2 k+l - v 2 k = -2(r k P k +x k Qk) - (r 2 k +x 2 k )^^ 



(1) 
(2) 
(3) 



where k = 0,---,N—l enumerates buses of the feeder, 
sequentially connected in a line, and Pk,Qk stand for real 
and reactive power flowing from bus k to bus k + 1. Vk is 
voltage, while pk and q k describe the overall consumption 
(negative/positive) of real and reactive powers at the bus 
k. The values of rk and Xk represent the resistance and 
inductance of the line element connecting k and k+l buses. 
Finite feeder line with the tap changing transformer in the 
beginning can be modeled by the boundary conditions in the 
rescaled p.u. variables: 



1. 



TV+l 



Qn+i = 0. 



(4) 



These boundary conditions close the system of equations 
( 1|2|3 1 and allow full reconstruction of the power and voltage 
profile given the consumption/injection, pk,q k - 

When the line is long and the number of consumers 
is large, N 3> 1, one is looking to simplify the system 
of equations and thus change to the continuous, N — » °°, 
limit. Assuming that all the power lines in the system are 
of the same grade, we set rk/xk to a constant. Then rk = 
rlk/L,Xk = xlk/L, where L is the total length of the feeder 
line. To develop the standard homogenization approach we 
assume that Pk,Qk,Vk are all decomposable into a sum of 
two components, one large but varying smoothly, and another 
small by amplitude but varying fast: F k = F(z) + F(L k )/N, 
where z = L k /L, and Lk = %!j=o ^ anc * both F and F are 
0(1). The homogenization approach suggests that, given 
homogeneous consumption along the feeder, p k and q k are 
small while p(z) — p k L/lk and q(z) — q k L/lk are 0(1). 
Relating differences to derivatives, F^+i — F k « F'(£)l k /L, 
and then applying it to all the relative differences in the Dis- 
tflow equations one arrives at the following set of Ordinary 
Differential Equations (ODEs) that will be called, naturally, 
DistFlow ODEs 




, P l +Q l 



rP+xQ 



(5) 



The DistFlow ODEs need to be solved under the following 
mixed (as defined on both ends of the feeder) conditions 



v(0) = l, P(L) = Q(L) = 0. 



(6) 



III. Inverter based control of reactive power 

Many nodes of existing feeders are capable to control their 
reactive consumption/production. This freedom is particu- 
larly handy in the modern PV systems which contain invert- 
ers [1]. Not only these inverters are capable of controlling the 
reactive power (obviously within their natural limits, usually 
related to the maximum level of real power the inverter is 
calibrated for) but they are also capable of doing so (at least 
in principle) in a smart and almost instantaneous manner, 
e.g. responding to some other measurements available for 
the feedback. A number of inverter-based control schemes 
responding to local voltage were recently discussed [1], [16]. 
Here, we will limit ourselves to analysis of the following 
two simple control schemes, which should all be contrasted 
against the "do nothing" case, when reactive power is in- 
jected by inverter in a natural proportion to the respective real 
power injection. Our first model (and very simple) control 
case is the one of the "zero power factor", when the reactive 
power flow from the node in the feeder is simply set to zero, 
q = 0. Another, and a bit more sophisticated, model analyzed 
here is similar to the EPRI example from [16], also discussed 
in [13], [14], [15]: 

q(v) = q (l - - ■ — 7 % ,^^1 , (7) 



l+exp(-4(v-l)/5) / 

where qo and 8 are positive constants, qo is related to control 
capacity of the inverter and 8 is the parameter expressing 
tolerance in the level of voltage variations. In this scheme 
with local feedback one assumes that q(z) responds to 
voltage measured at the same location, v(z). 

IV. Re-scaling and Simplification of the 
Boundary Value ODE Problem 

Assuming that p = const, and changing from the flow 
densities P, Q, voltage, v, and position along the feeder, z, 
to the new variables 



u(s) 



r P(z) 



\p\ v(L) 
v(L)' 




one arrives at the following rescaled version of the DistFlow- 
ODEs @ 




signO) - 2^3- 

p+Bx 
v 



(9) 



where A = q/\p\,B = x/r and the boundary conditions (j6j) 
turn into 



1>(0) = 1, p(0)=x(0)=0. 



(10) 



The rescaled formulation is advantageous as when q = 
const the resulting system of Eqs. ( 9|10 l defines a Cauchy 
problem, with the conditions ( 10 1 stated only on one end of 
the s, interval, in a pleasing contrast with the mixed original 
formulation of (5]6i. We solve the Cauchy problem ( 9|10 l 



forward in the rescaled time, s : — > s*, thus arriving at 
p(i*),x(s*),'V)(s*), and then one recomputes L, as well as 
P(0),g(0) and v(L), according to 

v(L) 



P(0) = 



pOOvIpIA 



»(*♦)■ 



e(0 ) = ^/p.(il) 



To summarize, Eqs. ( 9|10|1 1 1 allow efficient computations 
of the original mixed problem ( |5|6| l for different values of 
the feeder length L by simply scanning (increasing) s*. 

Note that the trick of reducing the original mixed problem 
Eqs. ( |9|10|11| > to single swap Cauchy computations gener- 
alizes to the special cases of the zero power factor control, 
q = 0. However, the technique does not apply to the case of a 
voltage dependent reactive control. Therefore, in discussing 
results of our simulations in the next Section we relay on 
the efficient Cauchy technique, but also (and when the latter 
does not apply) on resolving the boundary value problem 
directly. 

V. Case Studies of the DistFlow ODEs and 
Discussions 

This Section is devoted to analysis of the DistFlow ODEs 
in a number of interesting regimes. Utilizing the rescaling 
arguments above, one sets r—x—l in all of the simulations. 
We consider only lines constructed of statistically similar 
customers. Therefore, and in view of the homogenization 
and scaling arguments above, we will set the level of real 
power to p = qpl, depending on if the feeder is in the 
consumption or production regime. When no extra voltage 
control applies, this constant power consumption/production 
also translates into a constant consumption/production of 
the reactive power. We will set the base case of constant, 
q — p/2 reactive consumption/production and than compare 
it with two simple control schemes (assumed realizable by 
distributed inverters): zero power-factor control, q = 0, and 
feedback control on voltage described by Eq. Q. 

To set up the stage we, first, analyze the standard case, 
where both p and q are constant which are negative along 
the feeder, thus describing the feeder build of customers with 
constant (voltage independent) consumption. To generate 
data shown in Fig. ([TJ we fix p = — 1 and q = —0.5 and then 
scan different values of L. The result, shown on the top left 
sub-Figure of Fig. ([T]), is a standard in power engineering 
nose-curve illustrating the fact that the amount of power 
drawn from the system cannot exceed the threshold depen- 
dent on the system characteristics. Figure on the top right 
of Fig. ([T]) shows trade-off between the power-utilization of 
the feeder, measured in the ratio of the power injected at the 
head of the line, P(0) to the power consumed along the line, 
— p*L = L, and the voltage drop at the end of the line, v(L). 
Two bottom sub-figures in Fig. ([1} show profiles of voltage 
and power flows for lines with parameters correspondent to 
green and red markers in the top sub-figures. Markers select 
two representative stable solutions. (Here stability is defined 
with respect to the standard exciter dynamics, responding to 



p=-1, q=-0.5 




Position along the line, z g [0, L] — 



0.2 0.4 0.6 

Position along the line, ze [0, L] 



1 

0.8 
0.6 
0.4 
0.2 



1.1 
1 

i 

0.9 
• 0.8 
0.7 



p=-1,q=0 




0.2 0.4 0.6 0.; 
L 









(b) 



0.2 0.4 0.61 

Position along the line, z e [0, L] ^ 




0.2 0.4 0.6 

Position along the line, z e [0, L] 



Fig. 1. The case of uniform distributed consumption of real and reactive 
powers. 



Fig. 2. The case of uniform consumption and zero-power factor reactive 
control. 



the change in voltage. Then, the upper portion of the nose 
curve in Fig. where dv(L)/dP(L) > 0, is stable, while its 
lower portion, were dv(L)/dP(L) < 0, is unstable. We note 
that this definition of stability is a proxy to actual stability, 
which should be analyzed in the future in reference to actual 
dynamical equations, which are in fact different for different 
types of consumers/producers.) 

The bare case of distributed consumption, shown in 
Fig. (jlj, should be compared with Fig. (|2jl and Fig. d3J. We 
use the same type of diagnostics for creating the complex 
Figures as before, however we apply it now, in Fig. |2]) 
and Fig. (B}, to the case with the zero power-factor control 
(q = 0) and the reactive control of Eq. |7]i with feedback 
on voltage. Changing from the bare case to the zero power 
factor control does improve characteristics of the system by 
shifting the nose curve to the right and thus reducing voltage 
decrease for the feeder of the fixed length. The improvement 
continues with further transition from the zero power factor 
case of Fig. (|2]l to the case of reactive feedback control 
of Fig. p). These effects are quantitative and intuitive, 
and in fact not surprising as also following from an even 
simpler model with the distributed feeder replaced by a single 
clamped/aggregated load. 

Next we turn to analysis of a feeder with distributed 
generation, which may be though of as a homogenized 
proxy for a feeder rich on Photo-Voltaic systems generating 
enough of power during a sunny day to compensate for 
local consumption and injecting the access of power into 
the system. As before we showed in Fig. ( 4|5|6 1 the bare 
case (of comparable injection of real and reactive powers), 
the case of the zero power-factor (q = 0) control, and the 
case of reactive control with feedback to voltage respectively. 
We use the same diagnostics as above, showing voltage 
profile and power utilization as functions of the varying 
length of the feeder in the top sub-figures, and voltage and 
power flow profiles analyzed as functions of the position 
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Fig. 3. The case of uniform consumption and reactive control with feedback 
loop on voltage. 



along the feeder, for two (marked by green and red markers) 
representative cases. 

We observe that behavior in the Figures with distributed 
generation is qualitatively different from the one seen above 
in the case of distributed consumption. Most interesting new 
phenomenon, the one which was not seen in the case of 
distributed consumption, and which would not be seen in 
the naive model replacing the whole generating feeder by 
a single effective node, is associated with the emergence 
of multiple stable solutions. Red and green markers on the 
top sub-Figures of Fig. Q show solutions from the top (in 
terms of voltage) and the second top stable branches. Even 
though the lower voltage solution may seem unreachable 
as showing dangerously low voltage, the bare co-existence 
(for a range of L) of multiple stable solutions itself is 
extremely dangerous for the distribution system operations. 



Indeed, fault in the transmission system, resulting in an 
abrupt degradation of voltage to zero at the head of the line, 
may force the system to equilibrate in the post-fault regime 
to an undesirable low voltage branch. This scenario would 
be similar to the one associated with the phenomenon of 
the Fault Induced Delayed Voltage Recovery (FIDVR) [17], 
[18] which was observed in feeders rich on inductive motors 
(e.g. many air-condition loads with small inertia) [19]. The 
phenomenon of FIDVR can be interpreted as the one coming 
from a complicated hysteretic p(v) and q(v) characteristics 
of the individual inductive loads [20]. Similarly, we should 
worry about a hysteretic behavior in the case of a feeder 
with distributed generation. Notice, however that here the 
phenomenon does not originate from individual loads, but it 
rather emerges as a result of a complex spatial distribution 
of voltage and power flows along the feeder. 

We analyze details of voltage and power flow profiles 
for two representative feeders, selected in the top portion 
of Fig. Q by green and red markers. For the top branch 
(marked green) the profiles are standard: voltage raises 
monotonically, both real and reactive power flows are in- 
jected on the top of the line and then decrease monotonically 
(by absolute values) as we move along the line. However, the 
profiles are distinctly different in the second regime (marked 
by red markers in the top sub-Figures). Here, the direction 
of real flow reverses in the middle of the feeder - instead 
of flowing towards the head of the feeder (which would be 
the naively anticipated case with the export of power from 
this part of the feeder) the flow changes to local import, with 
the power flowing towards the end of the feeder. Note that 
this oscillatory behavior may be even more complex (with 
multiple reversals of the flow direction along the feeder) for 
other lower voltage branches. 

The two cases of control, applied to this bare case with 
distributed generators are shown in Fig. <|5j and Fig. (|6j. The 
changes in these cases (when compared with the bare case) 
are visible, however not alternating the principle observation 
made above: multiple stable solutions may still be realized 
in the feeder. We can understand this observation as follows. 
The main design suggestion for these schemes were based on 
considerations of voltage control. In this regards, the controls 
do achieve the anticipated goal as they do reduce increase 
of voltage in the top, base solution. However, the controls 
do not change the potentially dangerous co-existence of the 
reasonable (in terms of its voltage and power flow profiles) 
solution with its complex and very undesirable low-voltage 
sibling(s). 

One may question validity of modeling the distributed 
generation by generation with output independent of voltage 
in the low-voltage regime. We have tested a more realistic 
case (regularized at low voltage) as well (not shown in 
Figures) and observe that even if p decreases to zero with v in 
some reasonable manner, the main conclusion drawn above 
within the constant generation model stays - multiple stable 
solutions may emerge and standard reactive control, taking 
care of the voltage problems in the top branch, does not 
remove the operationally dangerous multiplicity/degeneracy 



p=1 , q=0.5 




0.5 1 1.5 2 1 0.5 1 1.5 2 



Position along the line, z e [0, L] ^ Position along the line, z e [0, L] 
o. 

Fig. 4. The case of uniformly distributed and comparable generation of 
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Fig. 5. The case of uniformly distributed generation of real power and 
zero reactive power (zero power factor) control 



of the power flow solutions. 

VI. Conclusions and Path Forward 

The key findings of this paper can be summarized as 
follows: 

• We have presented a novel approach to analysis of a 
long feeder with statistically homogeneous distribution 
of loads. When the number of consumers/buses is large, 
equations describing power flows over the feeder can 
be reduced to a system of ordinary differential equa- 
tions. The approach is computationally advantageous 
as containing fewer number of parameters than in the 
original discrete formulation. The technique also allows 
to account for distributed reactive and active power 
controls of individual loads. 

• In the (smart grid) future we expect some feeders to 
perform in the new regime of distributed generation. 
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Fig. 6. The case of uniformly distributed generation of real power and 
reactive control with feedback on voltage. 



We show that this regime exhibits a number of features 
distinctly different from the standard feeder of today 
containing only distributed consumption. We observe 
that in the distributed generation regime the power 
flow may reverse its direction along the feeder multiple 
number of times. Moreover, under certain conditions bi- 
and multi-stability can be observed, thus resulting in a 
hard to control, and dangerous, hysteretic behavior. 
• The distributed control of reactive power that is possible 
with programmable inverters provides principal capabil- 
ity of modifying the solution manifold. This option can 
be used to recover from an abrupt voltage fault as well 
as for trading off between losses and voltage deviations 
in the system. 

Note that the model(s) considered in the manuscript rep- 
resent only rough qualitative proxies for real power sys- 
tems with distributed generation. In spite of these models 
simplicity they show rich and nontrivial solutions, thus 
providing a short glimpse into the emerging complexity 
associated with the large-scale integration of distributed 
generation in modern power systems. More realistic and 
detailed models of power distribution systems should account 
for non-trivial voltage dependence of aggregated loads and 
generators on every bus. Accurate modeling is especially 
important for describing recovery from the transmission level 
induced faults in the feeder which was operating in the 
regime of distributed generation prior to the fault. It will 
also be important in the future to have a detailed dynamic 
analysis of the post-fault transients, in particular accounting 
for intermittent fluctuations of the distributed renewable 
sources. Finally, we anticipate that the technique introduced 
in this manuscript will be instrumental for addressing the 
fundamentally important (but still widely open) questions of 
formulating and designing optimal, centralized or distributed, 
control of reactive and active power flows in the power 
distribution systems. 
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